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The tidal gravitational field of galaxy clusters causes a coherent distortion of the images 
of background sources. Since the distribution of image distortions, namely the shear field, 
traces the local gravitational potential of a deflector, it can be used to reconstruct the two- 
dimensional mass distribution of clusters of galaxies. Moreover, the shear field can provide 
unique information on the redshift distribution of high-redshift galaxies. In this review we 



' summarize recently-developed parameter-free methods of cluster-mass reconstruction based 

on the shear field, and we apply a mass-reconstruction method to the cluster Abell 370 at 
£^ ! redshift 0.375. 

■ §1. Introduction 

o : 

Clusters of galaxies are the most massive gravitational bound systems in the 
universe and therefore contain crucial information on the origin of structure forma- 
tion and on cosmology. ^ " 3 ) For example, statistical studies of the cluster mass 
distribution and its evolution give us information on the distribution of dark matter 
and on the cosmological density parameter Qq. 

■ To this time, various methods have been used to determine the mass distributions 
in clusters of galaxies, such as: a dynamical method in which the velocity dispersion 
of member galaxies obtained from optical observations is used on the basis of the 

-£3 . virial theorem, 4 ) and a method based on diffuse X-ray emission from the hot (~ 

several keV) intra-cluster medium (ICM) that traces the gravitational potential of 
the cluster. 5 ) These two methods, however, are based on strong assumptions which 
are not fully justified. In applying the dynamical method, we must require that 
clusters are in virial equilibrium. However, this requirement is not always satisfied, 
because the typical dynamical time scale of a galaxy cluster is not much shorter 
than the Hubble time Hq , at least for a high-i?o universe, and the existence of 
substructures observed in some clusters indicate that they are not relaxed systems. 
Moreover, projection effects of substructures and the anisotropy of galaxy orbits in 
clusters may cause substantial errors in mass determinations by this method. On the 
other hand, X-ray analysis requires the assumption that the ICM is in hydrostatic 
equilibrium with the gravitational potential and also requires spherical symmetry 
for the mass profile of the cluster. These two methods are thus very sensitive to the 
physical and the dynamical state. 

Another approach to determine the mass distribution of clusters is to make use 
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of gravitational lensing, which can in principle be used to determine the projected 
mass distribution of the lens object directly, independent of the physical state and the 
nature of the matter content. The light rays from background sources are deflected 
by clusters of galaxies, and the strength of the deflection depends on the total mass 
and the mass distribution of the deflector. 

In the case of strong lensing, which occurs near the cluster core, the light-ray 
bending leads to the formation of arcs in cluster fields. If we assume spherical 
symmetry for the cluster mass distribution, the distance of an arc to the cluster 
center yields an estimate for the mass enclosed by a circle traced by the arc. An arc 
therefore places strong constraints on the cluster-mass distribution, especially in the 
central region of the cluster. 

On the other hand, weakly-distorted images of background galaxies — so called 
arclets — can be observed to much larger angular separations from the cluster cen- 
ter, and thus they provides valuable information concerning cluster-mass distribu- 
tions out to large radii. 6 ) Tyson, Valdes and Wenk 7 ) were the first to detect the 
weakly-distorted coherent images of galaxies behind two rich clusters. Kochanek 8 ) 
and Miralda-Escude 9 ) attempted to fit parameterized mass profiles of the clusters 
to the observed distortion fields. Kneib et al. 10 ) placed strong constraints on the 
bimodal mass profile for the cluster Abell 370 using one giant luminous arc (here- 
after referred to as GLA) and multiple images. Moreover, it has been shown that 
arc statistics are very sensitive to cosmological parameters, such as the cosmological 
constant A. n )' 12 ) Bonnet, Mellier and Fort 13 ) measured the coherent gravitational 
shear induced by the cluster CI 0024+16 out to 3 Mpc, where /150 is the Hub- 
ble parameter in units of 50 km s _1 Mpc -1 . It was first found by Kaiser and 
Squires 14 ) that the distortion field can be used for a parameter-free reconstruction 
of the two-dimensional mass density of a cluster in the weak-lensing regime. That 
is, the two-dimensional mass density of a deflector can be expressed as a convolution 
integral of the gravitational shear — which can be obtained from the image distor- 
tions of background galaxies in the weak-lensing limit — with a known kernel. The 
reconstruction formula derived by Kaiser and Squires was then generalized to include 
the strong-lensing regime 15 ) " 17 ) and modified to remove the artificial boundary ef- 
fects, 15 )> 18 ) " 21 ) which arise in the resulting mass map obtained with the original 
Kaiser and Squires mass-reconstruction technique. These mass-reconstruction meth- 
ods have been used to obtain the mass distributions of clusters in recent years (e.g., 
Refs. 22) ~ 25). 

In the present paper, we summarize the recent progress on the theoretical and 
the observational front in this field. There have been several reviews of this sub- 
ject. 26 ) ~ 28 ) We shall try to be self-contained and pedagogical so that the reader can 
perform the cluster-mass reconstruction by reading this paper. The remainder of 
the paper is organized as follows: In § 2 we present basic equations and concepts 
of gravitational lensing. In § 3 we introduce the definition of the image shape, and 
describe how the local observables are related to the local properties of the lensing 
clusters. In Section § we review some mass-reconstruction schemes that are based 
solely on image distortions. As mentioned earlier, applications of the reconstruction 
formula derived by Kaiser and Squires to real data involve several difficulties. Here 
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we describe the difficulties which one encounters in applying the original reconstruc- 
tion method to real clusters and how such difficulties can be overcome. Section 5 
summarizes the observational studies of galaxy clusters based on the shear analysis. 
In § 6 we apply a mass-reconstruction method to the cluster Abell 370 to demon- 
strate the feasibility of the shear-based analysis. There we employ a new method to 
break the mass degeneracy associated with the mass reconstructions based solely on 
image shapes. Finally, we conclude with a summary in § 7. 



§2. Gravitational Lensing 



In this section we summarize the basic equations and concepts involved in gravi- 
tational lensing which we shall need later. For general treatments, see Refs. 29) -31). 

2.1. Basic relations 

The gravitational field of a deflector changes the source position (5 to the appar- 
ent position 6 according to the lens equation. Let S{9) be the surface mass density 
of a deflector. The lens equation is written in terms of the two-dimensional effective 
lensing potential ip(9) in the form 

P = 9-V e if}{9). (2-1) 
Here if} (9) satisfies the two-dimensional Poison equation, 

A e iP(6) = 2k{9), (2-2) 
where k is the convergence defined by 

<¥) := (2 . 3) 

^cr 

with the critical surface-mass density 

Here D A , D s and D As are the angular-diameter distances from the observer to the 
deflector, from the observer to the source and from the deflector to the source, 
respectively. For a fixed lens redshift, we see that if>(9) depends on the source 
redshift through the distance ratio D<i s /D s in S CT . 

Since if} (9) satisfies the two-dimensional Poison equation (2-2), it can be written 
in terms of k as 

if>(9) = - f d 2 6' In \9- 6'\ k{9'). (2-5) 

7T J 

We then introduce the components of the shear, 

7i : = -^,22) 5 72:= "0,12, (2-6) 
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where an index i (i = 1,2) preceded by a comma denotes partial derivatives with 
respect to 9 l . It is useful to combine the shear components to define the complex 
shear, 

7 := 71 + i72 = |7| exp(2i^). (2-7) 
We further define the tangential shear 7t(#, 9') relative to the point 6' by 

7t (0; 9') := - (71 {9) cos [2ip(0, 9')] + 72 (0) sin [2ip(0, 9')] ) 



7 e 



-2»(v>+5) 



I7I cos 



2((f>-<p 



7T . 



(2-8) 



where denotes the real part of a complex number z, and ip(9, 9') is defined by 



f (9, 9') := arctan 



/1 - 



(2-9) 



From the definition, 7(6*) is a linear combination of the second derivatives of 4>{9), 
so that the complex shear can be expressed by the convolutional integral of k with 
a known kernel: 



7(0) 



'd\-d\ 



- I d 2 9'V{9-9')K{9'). 

7T J 



(2-10) 
(2-11) 



Here ft := d/d9 % (i = 1, 2) and V(9) is the complex kernel, 



V{9) :-- 



+ id 1 d 2 ln\9\ 



?2 — Q\ — 2i9\9 2 



(2-12) 



\e\* (9 1 -i9 2 ) 2 ' 

We note that there exists a global transformation that leaves 7(6*) invariant such that 

k{6) k{9) + k , (2-13) 

where kq is an arbitrary constant. We see in § 4 that this degeneracy leads to an 
ambiguity in the mass distribution predicted by the reconstruction technique based 
solely on the shear field. 

2.2. Magnification and image distortion 

The local properties of the lens equation are described by its Jacobian matrix: 



89 



-"0,12 
-Y>,12 1 - ^,22 



(2-14) 
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From Eqs. (2-2), (2-6) and (2-7), the Jacobian matrix A can be rewritten as 



A{0) 



1 



= (1-K 



and it has two eigenvalues 




cos(20) 
sin(2<£) 



sin(20) \ 
- cos(20) J ' 



At :=1-k±| 7 |. 



(2-15) 



(2-16) 



We see from Eq. (2-15) that the convergence k is the trace part, and the complex 
shear 7 is the trace-free part of the Jacobian matrix A(8). The convergence term 
alone in Eq. (2-15) causes an isotropic focusing of the light ray, while the shear 
term induces anisotropy in the lens mapping. The quantity \j\ = (jf + tI) 1 ^ 2 is the 
amplitude of the shear, and the phase <fi is the position angle of the eigenvector that 
belongs to the eigenvalue 7l_ = 1 — k — |— y | . 

The magnification of an image is given by the convergence and the shear as 
follows: 

1 (2-17) 



The closed curves defined by 




|deU| |(1-k) 2 -| 7 | 2 | 



= detA{9) 



(2-18) 



are called 'critical curves', on which 
the magnification diverges. The image 
plane is separated by critical curves; the 
regions where sign(det„4) = +1 and — 1 
are called even- and odd-parity regions, 
respectively. An intrinsically-circular 
source is transformed to an ellipse with 
axis ratio (< 1) of |yl_/yl+| for k < 1 
and for k > 1 and magnified 

by the factor fi = l/\A—A+\. The grav- 
itational distortion vanishes along the 
curve defined by n(9) = 1, which lies 
in the odd-parity region. 

In particular, for a non-singular 
circularly-symmetric lens, there are two 
circular critical curves: the inner critical 
curve defined by = A + = 1 — k + |7|, 
and the outer critical curve defined by 
= yl_ = 1 — k — \ Images close 
to the outer critical curve are elongated 
in the tangential direction, while images 




Fig. 1. Image-distortion field of intrinsically- 
circular sources for a non-singular circularly- 
symmetric lens. The dotted line represents 
the inner critical curve, where images are 
radially elongated. The dashed line repre- 
sents the outer critical curve, where images 
are tangentially elongated. Image distor- 
tions vanish along the curve k = 1 indi- 
cated by the solid line. The curve n = 1 
lies in the odd-parity region. 
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close to the inner critical curve are radially elongated (see Fig. 1). From these im- 
age properties, the outer and inner critical curves are called tangential and radial, 
respectively. 

§3. Observables vs lensing properties 

As stated in § 2.2, the tidal component of the gravitational field causes a co- 
herent distortion of the images of background sources. Hence, the galaxy images 
observed through a cluster can be used to infer the lensing properties of the cluster. 
For this purpose, we must first quantify the shape of a galaxy image in terms of 
observable quantities. For an elliptically-shaped source, its shape and orientation 
can be defined in terms of the axis ratio and the position angle of the major axis, 
respectively. In this case, an appropriate image ellipticity can also be defined by 
them. However, observations of galaxies reveal quite irregular shapes which can- 
not be accurately approximated by a simple ellipse, and this irregularity leads to 
serious noise in measuring the shear field. We must, therefore, take account of the 
irregularity in quantifying a galaxy shape. 

First, we define the center of a galaxy image by 

/«![!(«)] 

where 1(8) is the surface-brightness distribution of the source galaxy, and q(I) is a 
weight function, which we choose appropriately. For example, a possible choice for 
q(I) may be the Heaviside step function q(I) = H(I — 7 t h), where 7 t h is a limiting 
surface brightness, and H(x) is defined by 



*w H Z : t ■ (3-2) 




In this case 9 C is the center of the area enclosed by a threshold isophote 1(8) = 7 t h- 
For q(I) = /, 8 C is the center of light. We then introduce the tensor of second 
brightness moments for each image: 

^. Wim .^M. (3.3) 
Jd 2 8q[I(9)\ 

With this definition, we can quantify the shape of a galaxy image using the complex 
ellipticity e = e 1 + ie 2 : 19) ' 32) 

. = On ~ Q22 + 2iQ 12 = ,. 2i<p . . 

' Qii + Q 22 + 2(Q 11 Q 22 -Q2 2 )i/2 "I I • ^ ) 

If galaxy images have elliptical isophotes with axis ratio r (< 1), then the complex 
ellipticity is e = (1 — r)/(l + r)exp(2z^>), where the phase (f> G [0, it) is the position 
angle of the major axis. 
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Next, we consider the intrinsic (unlensed) surface-brightness distribution of the 
galaxy image I s (0), where the source position (3 is related to the image position 
9 through the lens equation (2-1). Liouville's theorem tells us that gravitational 
lensing conserves surface brightness, and therefore we can set 

h0) = I0). (3-5) 

We then define the tensor of second brightness moments Q si j of an unlensed source 
as 

" !<Pi3q\hW)} 

where /3 C is the angular position of the center of the unlensed source. If we assume 
that source images are smaller than the angular scale where the lensing properties 
change, then we can locally linearize the lens equation to obtain 

W-Pc)i = Aij(9 c )(e-9 c y. (3-7) 

From Eqs. (3-3), (3-6) and (3-7), we find that the tensors of second brightness 
moments of the source and image are related according to 

Q s = AQA T = AQA, (3-8) 

where A = A(9 C ) is the Jacobian matrix of the lens equation at angular position 9 C . 
If we further define the complex ellipticities e s of the unlensed source in the same 
way as in Eq. (3-4), the transformation between source and image ellipticity is given 
by 

e(e a) -I ^ + 9 ) / (I + 9* e s ) for \g\ < 1, 
etes ' 5j -\(l+ 5e *)/(4+5*) for \g\>l, 

where g{9) is the reduced shear defined by 

(j:=7/(l-«), (3-10) 

and the asterisk denotes complex conjugation. Since det„4 = (1 — n) 2 (1 — | <7| 2 ) , the 
condition \g\ < 1 (\g\ > 1) is equivalent to the condition det„4 > (det.4 < 0). 

Hence, the lensing properties are related to the observable quantity e only 
through the reduced shear g. In particular, the intrinsic ellipticity e s vanishes for 
circular sources, and thus we have 

t — i i — 

(3-11) 

In general, on the other hand, source galaxies have intrinsic ellipticities, so that only 
one image gives us no information on the local gravitational field. We can, however, 
extract the local lensing properties from the observed distortion field by assuming 
that source galaxies are randomly oriented. If this assumption holds, that is, if 
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e s = |e s | exp(2i(f)) has a random phase <p, then the probability distribution function 
Pe s (e s ) of the source ellipticity e s in general takes the form 

p es (e s )d 2 e s =p(\e s \)\e s \d\ts\^, (3-12) 

where the function p(|e s |) is normalized so that 1 = cZ| e s | |e s |p(|e s |). It has been 
shown by Seitz and Schneider 32 ) that the expectation value of the n-th moment e n 
is related to the reduced shear g through 



(e n ) es := Jd 2 e sPes e n (e s ;g) 

= J\\e s \\e s \p{\e s \)j^-e n {e s -g) 
J g n for \g\ < 1, 



~ } l/g* n for \g\ > 1. (3 ' 13) 

Note that the expectation value (e n ) €s does not depend on the source-ellipticity 
distribution p(|e s |); if we adopt a different definition of the image ellipticity (e.g., 
Ref. 16)) its expectation value depends on its distribution. 

In practice, however, we must replace the expectation value (e) £s with the average 
over a local ensemble of image ellipticities, e: 

(e) es (9)*e(9), (3-14) 

where 

Here A^gai is the number of galaxy images, and 6i is the angular position of the i-th 
galaxy (i = 1, 2, • • • , N) defined by Eq. (3-1), and W(9) is the weight function 

W(0) = ex P (-|^) (3-16) 

with smoothing scale A9. This scale must be small enough so that the lensing 
properties can be assumed constant over the effective smoothing disk of area ~ ttA9 2 
but large enough so that the smoothing disk contains a sufficient number of galaxies 
to suppress the random noise. 17 -*' 20 ) 

We can thus make use of the smoothed ellipticity field e(9) as a direct estimator 
for g(9) or l/g*(9). In the case of weak lensing (k ^ 1 and I7I ^ 1, i.e., |g| <^ 1), 
we have 

7 g e. (3T7) 



In this limit, the shear 7 is a direct observable. 
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§4. Cluster mass reconstruction by weak shear field 

In this section we review several cluster-inversion methods based on the shear 
analysis. The point of this study is that the two-dimensional mass distribution of a 
cluster can be directly obtained only from the observed shear field, and that we need 
not assume a priori a mass profile of a cluster. In § 4.1 we review the original mass- 
reconstruction method developed by Kaiser and Squires 14 ) and summarize several 
difficulties we encounter in a practical application of this method to real data. In § 
4.2-4 we describe how these difficulties can be overcome. In § 4.5 a method to infer 
the total mass within a circular aperture is described. 

4.1. Kaiser and Squires algorithm — linear inversion Formula — 

Since both k and 7 are linear combinations of the second derivatives of ip, using 
the relation in Fourier space, one can obtain an expression for k in terms of the 
complex shear 7; that is, the relation (2-11) can be inverted. 14 ) 

To see this, we express the convergence k{9) by its Fourier transform k(k) as 



K 0) = — / & k k $) ex P(^ • (4- 1 ) 



(27T) 

In Fourier space, Eq. (2-11) reads 

m = -k(k)V(k), (4-2) 

7T 

where k(k) and V(k) are the Fourier transforms of n(0) and T>{6), respectively, and 

they are defined in the same way as in Eq. (4-1). The Fourier transform of the 
complex kernel V(9) defined by (2-12) is 

V{k) = * kl ~ kl ^ 2iklk2 - (4-3) 
|fc| 2 

From Eq. (4-3), we find that T>T> = it 2 , implying 

fT 1 = 71-- 2 'V* ". (4-4) 
Inserting Eq. (4-4) into Eq. (4-2) yields 

k(k) = -j(k)V*(k). (4-5) 

7T 

This relation can be converted into the relation in the real 0-space as follows: 
K 0) -kq = - I d 2 e'V*(9- 9')j(9') 

7T J 

= - [ d 2 9'^[V*{6- 9')j(9')]. (4-6) 

7T J 

Here the constant kq is inserted into the right-hand side of Eq. (4-6) because adding 
a uniform mass sheet does not change the shear (see Eq. (2-11)). In terms of the 
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tangential shear j t defined by Eq. (2-8), the inversion equation takes the following 
form: 

K (0)- Ko = 1 -fd^^m. (4-7) 

V ' W \e-o>\ 2 

Thus we obtain an expression for the convergence n(9) in terms of the complex shear 
7(#) on the whole image plane. 

To apply the original Kaiser and Squires algorithm to real data, the following 
assumptions must be made: 

(a) The cluster is linear in the sense that k <C 1 and |7| ^ 1 everywhere. 

(b) Observational data are available over the entire space. This assumption may be 
dropped if the cluster is weak and small compared to the data field. 

(c) All background galaxies have the same effective redshift; i.e., all galaxies have 
approximately the same distance ratio D& s /D s . 

If these assumptions hold, we can perform the direct mass reconstruction from 
the observed distortion field e(9). In the case of weak lensing, the shear 7 is directly 
observable: 

7(0) « 10). (4-8) 

In practice, one has to replace the integral in Eq. (4-6) by a sum over a grid 6ij on 
which the shear is estimated. Using E{6) = E ct k(6), the surface-mass density E(9) 
of the lensing cluster is given by 

E(9) -Eq^E cx - I d 2 6' ^.[V* - 6') e0)} (4-9) 

TT J 

^E cr -J2 ^[D*0- %)e(%)], (4-10) 



7T 

h3 



where Eq is a constant surface-mass density which comes from k,q in Eq. (4-6), and a 
is the separation of the grid points. In this way, we can reconstruct the cluster-mass 
distribution from the weak-shear field up to an additive constant Eq. 

The crucial result here is that the projected- mass distribution of a cluster can be 
determined only from the observed shear field, that is, only from information about 
the shapes of galaxy images. However, we must require rather strong assumptions 
in deriving Eq. (4-10). Here we enumerate several difficulties which we encounter in 
performing the cluster-mass reconstruction with the Kaiser and Squires algorithm. 

(1) Application to non-linear clusters: Since the shear is observable only in the 
weak-lensing regime, the inversion algorithm described above cannot be applied to 
non-linear clusters. 

(2) Finite size of the data field: The integral in the inversion formula extends over 
the entire two-dimensional space, while real data are available only in a finite region 
restricted by the CCD area. Because of the lack of information outside the data field, 
we are forced to set 7 = there, which is equivalent for circularly-symmetric clusters 
to vanishing total mass within the data field. This sharp cut-off of the integration 
leads to an unphysical negative mass density near the boundary. We thus suffer from 
artificial boundary effects in mass reconstructions owing to the sharp cut-off. 
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(3) Redshift distribution of source galaxies: Cluster-mass reconstructions depend on 
the assumed redshift of the background galaxies through the distance ratio D& s /D s 
in U cr . Since the redshifts of background galaxies are unknown, there is a scaling 
ambiguity in the obtained mass distribution. Moreover, if background galaxies are 
distributed in redshift, the assumption that the distance ratio D^ s /D s is the same 
for all background galaxies does not hold, especially for high-redshift clusters. 

(4) Degeneracy of the solution for the cluster-mass inversion: We have seen that 
there exists a global transformation that leaves observable j(0) unchanged in the 
weak-lensing limit such that n(0) — ► k{6) + ko, with an arbitrary constant /-co; that 
is, the surface-mass density is determined only up to an additive constant from 
the observed shear field. We see in the next subsection that there exists a general 
global transformation that leaves the observable unchanged, which we encounter in 
all mass-reconstruction schemes based solely on image shapes. 

4.2. Generalization of Kaiser and Squires algorithm for critical clusters 

We have seen that the original Kaiser and Squires algorithm suffers from several 
difficulties. In this subsection we concentrate on the difficulties which we face in 
applying the Kaiser and Squires algorithm to non-linear clusters, where the shear 7 
is no longer a direct observable. 

Seitz and Schneider 17 ) generalized the Kaiser and Squires algorithm to include 
critical clusters, i.e., clusters which can produce critical curves. This technique, to 
be discussed in this subsection, is based on Eq. (4-6), and we need the assumptions 
(b) and (c) discussed in the preceding subsection. 

In order to extend the Kaiser and Squires method to the non-linear case, we 
express the inversion formula (4-6) with the reduced shear g. Using 7 = g(l — k) in 
Eq. (4-6), we have the integral equation 

K 0) = - [ d 2 9'V*(9- 0') g(0')[l - k(0% (4-11) 

7T J 

Here we have dropped a constant kq in Eq. (4T1) because the shear 7 is not an 
observable quantity in general. This equation can be formally expressed in an infinite 
power series as 

K {0)= g - Go g + go g Q g 

00 

= (4.12) 

n=l 

where g is the integral operator 

g(0,9') := - f d 2 6' V*(0-0')g{0') x, (4-13) 

IT J 

and g(9,9') acts on a function of 9'. The Kaiser and Squires method corresponds 
to the first-order approximation to this power series in the weak-lensing limit. 

For non-critical clusters, i.e., detA(0) > for all 0, the reduced shear is directly 
observable: g = (e) €s ~ e. For critical clusters, however, the relation between ellip- 
ticity and lensing property depends on the parity of the image (see Eq. (3-9)), and 
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we cannot determine the parity locally. To take account of the parity distinction, we 
write the shear 7 in the form 



7 = 3(1- «) 
= H(l - \g\) (1 - «) e + - 1) (1 - «) 4- 



(4-14) 



Seitz and Schneider 17 ) proposed an iterative procedure to solve the non-linear inver- 
sion equation (4-11): From Eqs. (4-11) and (4-14), we have 

= - f d 2 9'H(l - \g {n) (9')\) [1 - K (n \9')] 3? \v*{9- 6') e{9') 

7T J L J 



+ 



7T 



d 2 6' H(\gW(6')\ - 1) [1 - K (n) {9')}^ 



v*0-e') 



e*{9') 



(4-15) 



where gW = 7 ( n )/(l - k^) is the reduced shear in the n-th step of the iteration, 
and 7*™) is calculated by 



7 w = - f d 2 e'v{e- e') ^ n \e') 

vr J 



(4-16) 



This iteration is performed by starting with (6) = 7^°^ (6) = for all 6. Here we 
note that the first term on the right-hand side of Eq. (4-15) is the contribution to 
k{6) from the even-parity region, and the second term from the odd-parity region. 
In practice, of course, we use a discretized version of Eqs. (4-15) and (4-16). 

In general, the observable quantity is not the shear 7 but the reduced shear g 
(or l/g*). We see that the reduced shear g is invariant under the transformation 

k0) -► \k(9) + (1 - A) , 7(0) -► A 7 (^) (4-17) 

with an arbitrary scalar constant A 7^ 0; 16 - ) this transformation is equivalent to 
scaling the Jacobian matrix A{6) with A: 

' 1 - k - 71 -72 \ 
-72 1 - K + 71 I ' 



A{9) \A(6) = A 



(4-18) 



Thus we can determine k(9) only up to a constant. This is because we use the 
information only regarding the image shapes. We see that this transformation leaves 
the critical curves detA(9) = invariant. This indicates that we cannot determine 
the constant A even if GLAs, which are tracers of critical curves, are observed. 
Further, the curve k(9) = 1, on which the gravitational distortions disappear, is left 
invariant under the transformation (4-17). We can, however, constrain the value of 
A by requiring that the surface-mass density is positive everywhere, which yields a 
lower bound on k(9). 

A possible method to break this degeneracy is to employ the information about 
magnification effects. That is, we employ the fact that the magnification fj, transforms 
as 



(4-19) 
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under the transformation (4-17). Since gravitational magnification changes the size 
of a galaxy image and the number density of background galaxies, we can in principle 
break the mass degeneracy by combining the shear analysis with such magnification 
effects. 32 )" 34 ) 

4.3. Finite field inversion of the cluster mass distribution 

Both the inversion methods described above are based on Eq. (4-6) which in- 
volves the convolution integral of the shear over the entire two-dimensional space. If 
one tries to apply the inversion formula (4-6) to real data which are limited by the 
finite size of the CCD frame, then artificial boundary effects cannot be avoided. 

In this subsection we review the finite-field inversion method for the cluster- 
mass distribution derived by Seitz and Schneider. 20 ) This inversion method has been 
introduced to optimize the observational data, that is, to minimize the statistical 
errors due to noise. We assume, for simplicity, that all background galaxies have 
about the same distance ratio D^ s /D s as in § 4.1 and 4.2. We shall discuss the 
redshift distribution of background galaxies in § 4.4. 

All finite-field methods start from the fact that the gradient of the convergence 
k is related to the first derivatives of the shear 7: 15 ) 

V d K(8) = ( ™ + 72 ' 2 ) = U{9). (4-20) 
\72,i -7i,2 J 

In the weak-lensing limit where the shear 7 is a direct observable, we can obtain 
the vector field U(9) from the observed distortion field. However, the shear 7 is not 
an observable in general. Inserting 7 = g(l — k) into Eq. (4-20), we obtain the 
relation 15 ) 

V e K(e) = -±-( 1 -^ - 92 )( 9l ' 1 + 92 A=u(6), (4-21) 
l-\g\ 2 \-92 1 + gi J \ 92,1 -91,2) 

where K{9) is the scalar field defined by 

K{9) := ln[l - k0)\. (4-22) 

For a non-critical cluster, we can construct the vector field u{9) in terms of the 
observable quantity g{9) = (e) ts (9) ~ e(9). Hence, the relation (4-21) can be used to 
reconstruct the surface-mass density from the shear field on a finite region U. Since 
u{9) in Eq. (4-21) is a gradient field, the differential equations (4-21) can be solved 
up to an additive constant by integrating u(9) from an arbitrary point #0 to the 
point 9 along an arbitrary smooth curve I: 



K{9) = K(9 ) + all- u(l). (4-23) 



However, the vector field u{9) is obtained from noisy data, so that it will contain 
a rotational component. Therefore, u{9) is not a gradient field in general; in other 
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words, mass reconstructions from line integrals of Eq. (4-20) or (4-21) depend on 
the choice of the integration path. 

Seitz and Schneider 20 ) have explicitly taken into account the noise component 
of the vector field u(9) obtained from the observed distortion field by decomposing 
u(9) into gradient and rotational components: 

u0) = V e K{9) + V e x s(9), (4-24) 

where 

*x •<*>=- (4-25) 

Here s(9) is a scalar field which is assumed to account for the noise part within the 
data field U. The decomposition (4-24) can be determined uniquely by requiring 
that the average of Vg x s(9) within U vanishes, which is reasonable if the rotational 
component Vg x s(9) is due to noise. This requirement can be satisfied if s(9) satisfies 
the condition 

s(9) = const on the boundary dU of U. (4-26) 

(This is not a necessary but a sufficient condition.) If the condition (4-26) is satisfied, 
the scalar field K{9) can be obtained from the vector field (4-24): 

K{9) -K= [ d 2 9' H(9', 9) ■ u{9') (4-27) 
Ju 

with 

H(9', 9) := -%C{0', 9). (4-28) 

Here K is a constant which represents the average of K(9) within the data field U, 
and C{9', 9) is the solution of the Neumann boundary problem 

A e ,£(6', 9) = 5 2 (9- 9')-- V9eU (4-29) 

n{9') ■ %C{9', ff) = W G dU, (4-30) 

where A is the area of the data field U, and n is the unit vector orthogonal to the 
boundary dU; C(9',9) is solved uniquely up to an additive constant, so that the 
kernel H(9',9) is determined uniquely for a given boundary dU. The solution for 
two special geometries, a circle and a rectangle, is derived in Ref. 20). The kernel 
H(9', 9) is called noise filtering because it filters out the rotational noise component. 
We note that we can replace (K,u) in Eq. (4-27) with (k, U), because k and U are 
related to each other through the same relation (VgK = U) as that relating K and u 
(VgK = u); that is, 

K (9)-R= f d 2 9' H(9',9) -U(9>), (4-31) 
Ju 

where R is the average of k{9) within the data field U. We see from Eqs. (4-27) and 
(4-23) that K(9) = ln[l — k(9)] can be obtained up to an additive constant K; i.e., 
k(9) is determined only up to the global transformation 

k(9) -► Xk{9) + (1 - A) (4-32) 
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with an arbitrary constant A ^ 0. 

Several authors have developed inversion techniques +that require data only 
on a finite field. 15 -*' 18 - ) " 21 ) Seitz and Schneider 20 - ) have shown that there exist an 
infinite number of finite-field formulae, which are mathematically equivalent but 
different in their dependence on the noise due to the discreteness of data or the 
intrinsic ellipticity of the background galaxies. They also analyzed quantitatively the 
power spectra of the error fields for the mass maps reconstructed with synthetic data 
using different inversion methods: the non-linear version of the Kaiser and Squires 
inversion equation (4-15), the finite- field inversion equation with the noise- filtering 
kernel and two other finite-field inversion equations. The main results obtained there 
are (i) the inversion formula (4-27) performs better than the other two finite-field 
formulae, (ii) the noise for the finite-inversion (4-27) is more uniform than that for 
the other inversions over the data field and (iii) concerning the small-scale error, 
the modified Kaiser and Squires inversion (4-15) works best of all the inversion 
tested there, while it suffers from artificial boundary effects. They concluded from 
their results that the finite-field inversion with the noise-filtering kernel is the most 
accurate of all the inversion formulae used in their analysis. Lombardi and Bertin 
35 ) have obtained expressions for the errors on the mass maps reconstructed with 
various inversion formulae and shown that a rotation-free kernel, such as the noise- 
filtering kernel, should be used in order to optimize the reconstruction procedure. 
The analytical results obtained by Lombardi and Bertin confirm the conclusion of 
Seitz and Schneider. 20 ) 

4.4. The General cluster mass reconstruction algorithm 

Up to this point we have only dealt with the case where all background galaxies 
have the same effective redshift, i.e., all background galaxies have about the same 
distance ratio D& s /D s . This approximation is fairly good if the cluster redshift is 
relatively low and almost all galaxies have high redshifts. For high-redshift clus- 
ters, however, the redshift distribution of the source galaxies has to be taken into 
account in the inversion procedure explicitly. We review in this subsection the mass- 
reconstruction scheme developed by Seitz and Schneider, 32 ) which includes the case 
where background galaxies are distributed in redshift. Here we fix the lens redshift 
z& and express lensing properties explicitly as functions of the source redshift z (e.g., 
k(0) = k(0,z)). 

If we define the relative lensing strength w{z) for a source with redshift z by 

w(z) := H(z - z d ) pfL (4-33) 

with 

U CIiQO := lim E cr (z), (4-34) 
then the convergence and the shear for a source at redshift z can be expressed as 

k(9, z) = H(z - z d )^\ = w(z)p^- = w(z) Koo (0) (4-35) 
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and 

7 (6,z) = w(z) loo (8), (4-36) 

respectively. Note that k(9,z) = j(6,z) = for a source with redshift z < z&. The 
function w(z) can be calculated for a given background cosmology; for an Einstein-de 
Sitter universe (i?o = 1), w(z) is given by 

w(z) = H{z - z A ) -=== — , (4-37) 

\/l + z—\ 

and £ C r,oo is calculated to be 

cH (l + z d ) 2 
cr '°°~ SttGVT+z^-I 

= 1.39 x 1Q 14 ( 1 + Zd ) _ h5(j Mq Mpc- 2 . (4-38) 

V 1 + z d — 1 

From Eqs. (2-17), (4-35) and (4-36), the magnification n(9, z) of an image at angular 
position 6 and redshift z is given by 

MM ' 



|deU(0,z)| 

(4-39) 



[i-w(z) Koo(e)] 2 -w 2 (z)\ loD (e)\2 

If detA(0, oo) = [1 - KooiO)] 2 - hoo(0)\ 2 > for all 9, then the cluster is non-critical 
for all source redshifts. 

As we have seen in § 3, only g or l/g* is accessible to the observable quantity e. 
For a fixed source redshift z, the expectation value of e n (9) is 

(see Eq. (3-13)), where 

g(9, z) = 7( ^j = W(z)7 - (g \. (4-41) 
1 — k(9,z) 1 — w(z) k 00 (8) 

We then consider the case where source galaxies are distributed in redshift. Let 
p z (z) dz be the probability that the source redshift is within dz around z; this is a 
simple approximation, because the gravitational magnification changes the observed 
number density of galaxies, so that the probability distribution of the observed source 
redshift depends on the magnification fi(9,z). 33 ' This simplification is justified if 
the overall magnification is found to be small, or if the dependence of the galaxy 
redshift distribution on the flux is weak. 
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by 



If p z (z) is assumed to be known, the expectation value of e n can be expressed 



/•OO r 

(e n )e s ,z ■= J o dz Pz (z) J d 2 e sPts e n (e s ,g(z)) 



det4(z)>o 



dzp z (z)g n (z) + 



det.4(«;)>o 



dwp w (w) 



det„40)<o 
w 



dzp z (z) 



g*(z) 



loo 



\loo\ 2n JdetA(w)<o 



= % 



1 - KoqW 

dwp w (w) 
1 



1 - K 00 W \ H 
W 



^n( K 00)7oo) "i" I I n Y n (Kqo , 7oO ) 

l7oor n 



(4-42) 



where p w (w) is the probability distribution function of w, which is given by p w (w) dw 
p z dz. Here the functions X n and Y n depend on and 7^ by 



-X«,( k ocm 7o 



Jo 



min(l,l/(/t 00 +|7 00 |)) fl 
+ 



/_ 



xdw p w (w) 



w 



1 - KooW 



max(l,Kcx)-|7cx)l) ' 

n 



Yn^Koo, 7o 



max(l,Koo - |7oo I) 

min(l,l/(ft +I700D) 



dwp w (w) 



1 - Knr-W 



w 



(4-43) 
(4-44) 



Now we have the local relation between the expectation value of e and the lensing 
properties, 



7oc(0) = {e)e s ,M 



1 



l7ooW| : 



:H(«ooW,7ooW) 



(4-45) 



In a practical application, the redshift-averaged ellipticity (e) es;2 (0) must be replaced 
by the observed ellipticity e(6) defined by Eq. (3-15). In the weak-lensing limit (i.e., 

Koo <C 1 and \joo\ ^ 1)> we have 



^n(Koo,7oo) ~ {W n ) ; y n (Koo,7oo) =0 

with 

: = / (iu;p w (u;)u; n , 
JO 

and therefore the shear is observable: 

7oo$ « e(0)/<«>>. 



(4-46) 
(4-47) 

(4-48) 



That is, in the weak-lensing limit, the situation is the same as in the case where all 
sources have the same redshift z such that w(z) = (w). 
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In general cluster-mass reconstructions, the finite-field inversion formula (4-31) 
can be used: 32 ) 



Koo(0) - Koo = / d 2 e'H(9',e) • c4(#) 

Ju 



with 



(700)1,1 + (7 00 1 2,2 
(700)2,1 - (700)1,2 



(4-49) 



(4-50) 



Here /too is the (unknown) average of Koo(O) within the data field IA. Inserting Eq. 
(4-45) in Eq. (4-49), we have the integral equation for Koo(9), which can be solved 
iteratively. Starting the iteration with k& 0) = 7^? '(9) = for all 9, we calculate 
Koo and 7^ for n > 1 by 



K (^)(9)-R 00 = / d 2 9'H(9',9)-U<£ +1 \9>), 
Ju 



An+l) = - 
too c 



An) 2 
700 



-1 



(4-51) 
(4-52) 



Here the vector field Uoo +l ^ 



(9) is calculated from Eq. (4-50) with ^ +1) {9). 
It can be seen from Eq. (4-51) that this iterative procedure still contains the 
undetermined constant R^. Hence, this inversion algorithm involves a global trans- 
formation that leaves the observable invariant. In addition, the dependence of the 
resulting mass distribution on the constant Rqo is highly non-linear, so that this trans- 
formation cannot be determined analytically. In the case of weak lensing (kqo <C 1 
and <C 1), we have seen that the mass reconstruction depends only on (w), so 
that the global invariance transformation is equivalent to adding a constant. If the 
cluster is non-critical for all sources (i.e., det^4(oo, 9) > for all 9), then Y n vanishes 
and X n depends only on Kqo, yielding a somewhat simple relation: 



<6«> es , z = 7 ™ *„(««,). 
On the other hand, an approximation for the function Xi(koo) is 



^i(Koo) 



(w) 



1 — K n 



(w 2 ) 
(to) 



(4-53) 



(4-54) 



which is quite accurate for non-linear but not very strong clusters. 32 ) In this case, 
the global invariance transformation becomes 



^00 

(9) - \ Koo (9) 



(i-A)M 

(w 2 ) 



(4-55) 



with an arbitrary constant A / 0. 
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4.5. Cluster mass estimates from shear fields 

In the previous subsections we have described the cluster-mass reconstruction 
methods based on the shear analysis. One of the goals of such studies is to obtain the 
total mass of a cluster within a given region. In estimating the cluster mass inside a 
given aperture, one method to obtain the aperture mass is to use the reconstructed 
two-dimensional mass distribution. However, the resulting mass map from the shear 
field will be noisy because of the discreteness of the galaxy images, intrinsic source 
ellipticities, etc. In addition, it is difficult to estimate an error for the local conver- 
gence since the errors of the reconstructed convergence at different points will be 
strongly correlated. In this subsection we describe a method to infer the projected 
cluster mass inside circular apertures from weak lensing without reconstructing the 
entire mass distribution. 

First, we define a polar coordinate system 92) centered on a point #0 such 

that 

<4 - 56) 

Then the convergence averaged within a circle of radius # around 6q is given by 
k(0) := -4t r d& & I dtp' k(9($', if/)). (4-57) 

TTV 2 JO J 

Using the two-dimensional version of Gauss's theorem, this quantity can be trans- 
formed into 

«(<?) = <K>(0) + <7t>(0;&), (4-58) 
where the angular brackets denote the average over a circle of radius e.g., 

<7t)(tf;0o) :=/^7t(0>,^);0o) (4-59) 

with the tangential shear jt(9; 9q) relative to the origin 9q of the coordinate system. 
On the other hand, can be expressed as 

= A I* d#tf(K)(&). (4-60) 
v Jo 



From Eqs. (4-58) and (4-60), we have 

dR{$) _ 
dim? " 



-2<7t>(0;0 o ). ( 4 - 61 ) 



Integrating Eq. (4-61) between radii $ and ■&' (> ■&), we obtain the so-called (- 

statistic, 22 ) 

m tf) := m - tf') = -^—^ ^ — < 7t )(tf"; 0o), (4-62) 
where is the average of the convergence within an annulus between 1? and 

*{*,&) := W) [ d»" $"{„){$"). (4-63) 
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In the weak-lensing limit (i.e., k <C 1 and | — y | <C 1), the (-statistic can be deter- 
mined uniquely from the shear field, because the invariance transformation, which 
corresponds to an additive constant to k, cancels out on the right-hand side of Eq. 
(4-62). 

Since has a non-negative value, ((6,6') yields a lower limit on 

that is, the quantity TT(D d 'd) 2 E cr ( > ('d, "&') yields a lower limit on the lensing mass 
inside a circular aperture of radius Thus we can infer the lensing mass inside 
a circular boundary from the data outside the boundary in the weak-lensing limit. 
The (-statistic is useful in estimating the total cluster mass within the data field. 

We note that the shape of an aperture need not be restricted to a circle. Aperture 
masses for arbitrary aperture shapes are dealt with in Ref. 36). 

§5. Observational Studies of Weak Lensing 

This section reviews the observational studies of galaxy clusters based on the 
weak-lensing analysis. In § 5.1 we summarize some observational results from weak 
lensing. In § 5.2 we discuss the prospects of weak-lensing analysis on clusters with 
the upcoming 8.3 m Subaru telescope. 

5.1. Observational results 

The non-parametric methods for cluster-mass reconstructions by means of weak 
lensing have been widely used in investigating cluster properties in recent years. The 
pioneering work in this field was carried out by Kaiser and Squires, 14 ) who derived 
an exact inversion equation for the two-dimensional mass density of the deflector in 
terms of the tidal component of the gravitational field, namely the shear which is 
observable in the weak-lensing regime. The direct mass-reconstruction method of 
Kaiser and Squires and its variants were then applied to real clusters. In practice, 
these mass-reconstruction methods were mainly used to obtain the morphology of 
the mass distribution in a cluster. On the other hand, the (-statistic was often used 
in the cluster-mass estimations, as mentioned in § 4.5. We summarize in Table I 
some results of the weak-lensing analysis on clusters. 

Fahlman et al. 22 ) first applied the Kaiser and Squires method to the X-ray lumi- 
nous cluster MS 1224+20 at redshift 0.33. The shear field for this cluster was mea- 
sured from the observational data obtained with the 3.6 m Canada-France-Hawaii 
Telescope (CFHT), and the two-dimensional mass map was derived from the shear 
field. The location of the main peak of the resulting mass map is in good agreement 
with that of the optical and X-ray centroid. The mass estimates were performed 
using the (-statistic (4-62) within a circular aperture of radius 2.'76 (corresponding 
to a physical radius of 0.96 /i^ 1 Mpc for an Einstein-de Sitter universe) centered 
on the mass peak. They obtained a lower bound on the aperture mass interior to 
radius 0.96 h^Q Mpc to be ~ 7.0 x 10 14 K^qMq, which corresponds to a mass-to- 
light ratio of M/L ~ 400 h^Q in solar units. The mass estimate for this cluster from 
weak-lensing analysis is about three times larger than that from virial analysis. 

Smail et al. 37 ) analyzed three clusters which span a wide redshift range, z d €. 
[0.26,0.89]. Two of these clusters (CI 1455+22, z d = 0.26; CI 0016+16, z d = 0.55) 
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were selected for their high X-ray luminosities, and the other (CI 1603+43, z d = 0.89) 
was optically discovered and has a low X-ray luminosity. They detected clear lensing 
signals in the two lower-redshift clusters, while no significant lensing signal was 
detected in the highest-redshift cluster CI 1603+43. They concluded from their 
results that the absence of any lensing signal for their highest-redshift cluster is 
ascribed to the absence of a high-redshift population of background galaxies with 
/ < 25. On the other hand, Luppino and Kaiser 24 ) argued that the failure of Smail 
et al. to detect weak lensing in CI 1603+43 is simply due to the fact that this 
cluster is not massive enough to produce a detectable shear signal. The detection 
of significant shear signals in high-redshift clusters (MS 1054-03, z d = 0.83; 24 ) MS 
1137+66, z d = 0.783; 38 ) RXJ 1716+67, z d = 0.8 13 38 )) indicates that the great part 
of background galaxies are at redshifts of the order of unity, supporting the argument 
by Luppino and Kaiser. 

Seitz et al. 25 ) applied a direct mass-reconstruction method to the distant cluster 
CI 0939+47 at redshift z d = 0.41, observed with the Wide Field Planetary Cam- 
era 2 (WFPC2) on the Hubble Space Telescope (HST). Owing to the high redshift 
of this cluster and the small field-of-view of the WFPC2, the mass-reconstruction 
scheme described in § 4.4 — which takes into account the redshift distribution of 
the background galaxies and is based on the finite-field inversion equation — was 
used. Assuming a particular form of the parameterized distribution function of the 
source redshift, they reconstructed two-dimensional mass maps for several sets of the 
parameters. They found a strong correlation between the reconstructed mass map 
and the light distribution of the bright cluster galaxies. The main mass and light 
maximum agree well with a maximum in the X-ray image from the ROSAT/PSPC 
observation. 39 ) They also found a magnification effect on the observed number den- 
sity of background galaxies, by which they removed the mass degeneracy. 

Direct mass-reconstruction methods have also been applied to the QSO fields. 
From the CFHT observation, Fisher et al. 40 ) measured a weak gravitational shear 
induced by the cluster at a redshift of 0.355 in the field of the double QSO 0957+561, 
which is a multiply-lensed system with a well-studied time delay. They reconstructed 
a two-dimensional mass map in this field from the measured shear field in order to 
construct a detailed model for the lens system which consists of a primary lensing 
galaxy and the cluster containing this galaxy. The resulting mass distribution is 
consistent with a spherical cluster possessing a mass peak slightly offset from the 
position of the primary lensing galaxy. The obtained projected mass within a circular 
aperture of radius 1 Mpc centered on the mass peak is 3.9 ± 1.2 x 10 14 h^M Q . 

Several authors have compared cluster-mass estimates from weak lensing and 
those from other methods. Detailed quantitative comparisons of cluster-mass es- 
timates on scales ~ 1 Mpc show that mass estimates from weak-lensing analysis 
are consistent with those from X-ray analysis based on the hydrostatic equilibrium 

of the ICM with the gravitational potential (Abell 2218; 41 ) Abell 2390; 42 ) Abell 
2 163 43),44) ) _ 
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5.2. Observation with the Subaru telescope 

The weak-lensing analysis of galaxy clusters requires quite a large number den- 
sity (~ 50 arcmin -2 ) of galaxy images with sufficient detail and accuracy in order 
to obtain a reliable shear measurement. Hence, cluster-mass reconstructions based 
on the weak-shear field rely on the power of telescopes and instruments to be used. 
Moreover, for ground-based observations, the effect of atmospheric seeing — the cir- 
cularization of galaxy images — will reduce the lensing strength considerably since 
most of the background galaxies are quite small and faint. Hence, an excellent see- 
ing condition (sub-arcsec) as well as a high angular resolution is required for the 
ground-based weak-lensing analysis. 

On the other hand, the space-based deep HST/WFPC2 observation provides us 
with a large galaxy-number density without the seeing effect. Hence, it is a natural 
consequence that HST/WFPC2 observations have been frequently used in weak- 
lensing analysis in recent years. 25 )' 45 )' 46 ) However, in contrast to these advantages, 
the HST/WFPC2 observations have a fatal weakness from the theoretical point of 
view. As mentioned above, a mass reconstruction based only on the shear field suffers 
from mass degeneracy. If the cluster-mass distribution extends beyond the observed 
data field, we cannot uniquely determine the solution for the cluster-mass inversion 
from the shear field. In this case, we obtain only a lower limit on the cluster mass 
in the observed field. Moreover, since the dependence of mass reconstructions on 
the free parameter (see § 4.3) is highly non-linear for a high-redshift and super- 
critical cluster, even the resulting morphology of the cluster could be uncertain, 
depending on the choice of the parameter. In addition, the irregular shape of the 
WFPC2 frame will also make a mass reconstruction difficult to realize. Hence, both 
the requirements of a large galaxy-number density and a wide field-of-view are needed 
for a reliable mass reconstruction. 

The Subaru telescope is among such telescopes which ideally satisfy both the 
requirements. Subaru is the Japanese 8.3 m optical-infrared ground-based telescope 
at the summit of Mauna Kea, Hawaii. Among the observational instruments for 
Subaru, Suprime-Cam (Subaru Prime Focus Camera) is the most suitable for weak- 
lensing analysis. Suprime-Cam at the F/2.3 prime focus has a wide-field imaging, 
covering a 30' x 24' field-of-view with a 0."2/pixel resolution. 47 ) This corresponds 
to a physical scale of about 7.6 x 6.I/150 2 Mpc 2 at a moderate redshift of 0.2 for 
an Einstein-De Sitter universe. Weak-lensing analysis using the Subaru telescope 
with Suprime-Cam can probe the mass distribution on large scales (~ 10 Mpc) 
where the formation of gravitational structures is still in progress. Such studies can 
provide unique, invaluable information on the distribution of dark matter and Qq. 
For example, the super cluster MS0302+17 containing three clusters (CI 0303+1706, 
z d = 0.418; MS 0302+1659, z d = 0.426; MS 0302+1717, z d = 0.425 ) is of interest 
to investigate such problems. 4S ) 
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Table I. Summary of observational studies of clusters from weak lensing. 
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Notes: 

a) Redshift of the cluster. 

b) Radius of the aperture in units of Mpc: (*) area of the observed field. 

c) Projected cluster mass within the aperture of radius r in units of 1O 14 M0. 

d) Mass-to-light ratio of the cluster in solar units. 

e) Telescopes: i) Cerro Tololo Inter-American Observatory telescope; ii) Canada-France-Hawaii 
Telescope; iii) Hubble Space Telescope; iv) William Herschel Telescope; v) University of Hawaii 
telescope. 

f) References 
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§6. Application to the cluster Abell 370 

6.1. Background 

Abell 370 is a very rich, distant cluster of galaxies at a redshift of 0.375, domi- 
nated by two bright cD galaxies. To this time, one GLA, several multiple images and 
a number of arclets have been observed in this cluster. The GLA was discovered by 
Soucail et al. 52 ) and identified spectroscopically as the image of a background galaxy 
at a redshift of 0.724 lensed by Abell 370. 53 -* Since the discovery of the GLA, the 
cluster has been the subject of extensive lensing studies and observations. 10 )> 54 ) -58 ) 
Kneib et al. 10 ) investigated the mass distribution in Abell 370 on the basis of their ex- 
cellent ground-based CCD image. Assuming that the GLA consists of three merging 
images, they constructed a detailed mass model based on two elliptical components 
associated with the two bright cD galaxies by fitting the model to the observed 
GLA and multiple images. Smail et al. 55 ) discovered a radial-arc candidate in the 
HST/WFCl image. Using the Kneib et al. model, they predicted its redshift to be 
1.3 ± 0.2. Ota et al. 57 ) modeled the projected mass distribution in Abell 370 on the 
basis of the ASCA observation, assuming spherical mass profiles. They compared the 
radial profile of the projected mass inferred from X-ray analysis with that inferred 
from the Kneib et al. model based on strong lensing. The lensing mass interior 
to the GLA radius of 160/i^q 1 kpc is about three times larger than the X-ray mass. 
They attributed the mass discrepancy for this cluster to the projection effect of sub- 
structures. Recently, AbdelSalam et al. 56 ^ developed a non-parametric method to 
reconstruct the cluster-mass distribution based on the observational constraints by 
strong lensing. They applied their non-parametric method to Abell 370; they divided 
the projected cluster mass into squire pixels and searched for one of the solution for 
the mass distribution that follows galaxy light as closely as possible. From a visual 
inspection of the HST/WFCl image, 55 - ) AbdelSalam et al. interpreted the GLA 
as a five-image system in contrast to the Kneib et al. model. Using the GLA, the 
multiple images and the radial arc, they obtained mass maps for Abell 370, which 
reveal an obvious bimodal feature in the mass distribution with the two mass peaks 
close to the two cD galaxies. However, the resulting two mass maxima were slightly 
closer to each other than the two cD galaxies. In addition, the southern mass clump 
was found to be more massive than the northern one, though the northern cD galaxy 
is brighter than the southern one, consistent with the Kneib et al. model. The in- 
teresting result obtained there is the presence of an extra substructure close to the 
two mass peaks, which does not follow the light. Bezecourt et al. 58 ** improved the 
Kneib et al. model using the deep WFPC2 image. Detailed information concerning 
the GLA, multiple images and arclets placed strong constraints on the mass model 
consisting of the cluster- and galaxy-scale components. 

In this section we present a non-parametric reconstruction of the mass distribu- 
tion for Abell 370 based on weak-lensing analysis. In contrast to the mass recon- 
struction based on strong lensing as performed by Kneib et al., AbdelSalam et al. 
and Bezecourt et al., we do not use the information regarding the redshifts of the 
arcs: We make use of the information on the image shapes alone. Throughout this 
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section, we assume a matter-dominated Einstein-de Sitter universe with Qq 
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6.2. Observation and data reduction 

The data for Abell 370 have been 
retrieved from the HST archive. The 
cluster Abell 370 was observed in De- 
cember 1995 using the WFPC2 cam- 
era with the F675W filter (Rq75w) on 
the HST (ID: 6003, P.L: R. P. Saglia). 
The total exposure time is 5600 sec- 
onds. After STScI pipeline processing, 
the data were shifted and combined into 
the final frame to remove cosmic rays 
and hot pixels using the IRAF/STSDAS 
task CRREJ. We discarded the PC chip 
from our analysis because of its brighter 
isophotal limit, and thus the final frame 
consists of three WFPC chips. The side 
length of the data field is about 2.'5, cor- 
responding to O.93/150 Mpc at redshift 
Zd = 0.375. (1" on the sky represents 
6.2^ kpc.) 

To construct a catalog of faint ob- 
jects in this frame and measure the 
shape parameters (i.e., the center, the 
size and the complex ellipticity e defined 
by Eq. (3-3)) for each object, we used 
the SExtractor package. 59 - ) We selected 
all objects with isophotal areas larger 
than 12 pixels (0.12 arcsec 2 ) above a de- 
tection threshold of 2a /pixel, corresponding to /i675W = 24.5 mag arcsec -2 . The first 
and second brightness moments were computed for each object to determine the 
shape parameters. The faint and the bright magnitude limits were chosen so as to 
give reliable shape parameters. A catalog for arclet candidates was constructed with 
a total of 177 galaxies in a magnitude range i?675W £ (23, 25). 

As mentioned above, one GLA [z = 0.724) and one radial-arc candidate have 
been observed in the field of Abell 370. Through the detection process described 
above, the radial arc was identified as an object with i?675W = 24.0, and thus it 
was included in our catalog of faint galaxies with -R675W £ (23,25). On the other 
hand, the GLA was identified as a three-component image with i?675W = 20.7, 20.6 
and 20.5, so that this image system was excluded from our catalog of faint galaxies 
with i?675W £ (23,25). However, the GLA provides us with invaluable information 
on the mass distribution especially in the central region, and thus it will strongly 
constrain the cluster profile. Hence, we include this image system in our galaxy 
catalog. Thus the total number of the arclets contained in our galaxy catalog is 



x [arcmin] 

Fig. 2. Relative locations of the two cD galax- 
ies and the image systems contained in our 
catalog of 180 galaxy images. North is to 
the top, and East is to the left. The po- 
sitions of the two cD galaxies are marked 
with +. The positions of the giant lumi- 
nous arc and the radial arc are marked with 
A and □, respectively. The positions of the 
other arclets are marked with O. 
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A'gai = 180, corresponding to 38 galaxies arcmin" 2 . In Fig. 2 we show the relative 
locations of the two cD galaxies and the image systems contained in our catalog of 
180 galaxy images. 

6.3. Shear field 

In this subsection we describe the 
shear analysis of Abell 370. First, we 
divide the CCD-data field consisting of 
three WFC chips into two rectangular 
fields A and B; field A with a side 
length of about 2.'5 x 1.'25 consisting of 
the lower-left and the lower-right WFC 
chips, and field B with a side length of 
about 1/25 x 2/5 consisting of the lower- f 
right and upper-right WFC chips. Both j_, 
the cD galaxies are located at the left s 
edge of the lower-right WFC chip (see 
Fig. 2), and thus both the fields A and 
B contain them. For field A, the num- 
ber of the galaxies contained in our cat- 
alog is 142, corresponding to a galaxy- 
number density of 46 arcmin -2 . On the 
other hand, for field B, the number of 
galaxies contained in our catalog is 107, 
and the corresponding galaxy-number 
density is 35 arcmin -2 , which is much 
smaller than a typical galaxy-number 
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density of ~ 50 galaxies arcmin 2 in a 



weak-lensing analysis. This small num- 
ber density is ascribed to the small num- 
ber of faint galaxies detected in the 
upper-right WFC frame from our anal- 
ysis (38 galaxies in this frame; see Fig. 
2). 

Next, we introduce a rectangular 
grid with a constant grid separation of 
0/083 (Sl/igo 1 kpc) for each field: field A 
with 30 x 15 grid points, and field B with 
15 x 30 grid points. In order to reduce 
the noise due to the intrinsic ellipticities 
of the background galaxies, we calculate 
the local mean image ellipticity e(9ij) at 
each grid point 9ij from the image ellipticity e 



Fig. 3. The shear field obtained from 180 
galaxy images taken with the HST/WFPC2. 
The orientation and the length of a vector 
indicate the orientation <f> and the absolute 
value |ej of the locally-averaged image el- 



lipticity e = | e | e 2l< *, respectively. The case 
|e| = 1 (i.e., vanishing axis ratio) corre- 
sponds to a vector of length 0.'3. We adopt 
an optimal smoothing length AO depend- 
ing on the distance from the positions #n 
and #s of the northern and the southern 
cD galaxies, respectively; 6* N = (1'33, 1.'08) 
and S = (1/29,0/45). The smoothing 
length AO ranges from ~ O.'l (37 kpc) to 
~ 0.'3 (112 kpc), depending on the grid po- 
sition; 1' corresponds to 0.37 Mpc at the 
cluster redshift z<j = 0.375 for H = 50 km 
s" 1 Mpc -1 and an Einstein-de Sitter uni- 
verse. The positions of the two cD galaxies 
are marked with +. 

] k) of the fc-th galaxy at angular posi- 



tion 6k (k = 1, 2, ■ ■ ■ , iVgai; iVg a i = 180) (see Eqs. (3-15) and (3-16)). Here we employ 
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an optimal smoothing procedure; we adopt an optimal smoothing length A9(9ij) in 
Eq. (3-16) depending on the grid position 6%j such that 

A90ij) = min(0.12 |% - 0n| + A6 , 0.12 \9 {j -9 S \ + A9 ) arcmin, (6-1) 

with constant smoothing length A9q = 0/1, where #n and 6>s are the angular positions 
of the northern and the southern cD galaxies, respectively; #n = (1/33,1/08) and 
#s = (1/29, 0/45) in our coordinate system. Thus the smoothing length A9(9ij) at 
grid position 9ij ranges from ~ 0/1 (37/15Q 1 kpc) to ~ 0/3 (112h^ kpc), depending 
on the distance from the positions #n and 9s of the two cD galaxies. 

Figure 3 displays the resulting map of the locally-averaged image ellipticities 
e = |e| e 2 *^ obtained using 180 arclet candidates. The orientation of a vector indicates 
the direction <f> of the shear, and the length of a vector is proportional to the strength 
of the shear. A coherent shear pattern induced by the gravitational field of the cluster 
is clearly visible. 

6.4. Mass reconstruction 

Using the smoothed image-ellipticity field e(9ij) obtained from our galaxy cat- 
alog, we reconstruct the mass distribution of the cluster Abell 370. To take into 
account (1) the strong-lensing features in the cluster field, (2) the fairly high redshift 
of the cluster (z d = 0.375) and (3) the small field-of-view of the HST/WFPC2 field, 
we follow the generalized mass-reconstruction scheme 32 ) based on the finite-field 
inversion derived by Seitz and Schneider 20 ) (see § 4.3 and 4.4). 

Because of the irregular shape of the 
WFC field, we separately reconstruct 
the two-dimensional convergence fields 
,a(%) and Koo,b(%) on two rectan- 
gular fields A and B, respectively, as 
done in Ref. 25). Since a mass recon- 
struction with the finite-field formula 
(4-31) yields a constant Kqo, which rep- 
resents the average of Kco over the data 
field, we have two free constants in our 
reconstruction, Koo,A f° r field A and 
^oo,B for field B. One of the two con- 
stants can be used to join the two re- 
constructions, and the residual constant 
can be used to normalize the cluster 
mass. The procedure to determine the 
two constants and to join the two inde- 
pendent reconstructions is described in 
§ 6.5. 

Moreover, we must know the redshift distribution of the galaxy population in 
the mass reconstruction (see § 4.4). However, the redshift distribution of faint galax- 
ies with -R675W £ (23, 25) is uncertain, so that we adopt a parameterized redshift 




Fig. 4. The solid line indicates the relative 
lensing strength w(z) defined by Eq. (4-37) 
for the cluster redshift z<\ = 0.375. The 
dashed line indicates the assumed redshift 
distribution p z (z) of background galaxies 
defined by Eq. (6-2) for /3 = 1 and (z) = 
3z = 0.7. 
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distribution p z (z) of field galaxies of the form 25 )' 60 ) 

p * ) = T^ ffip (-<*/<)> < 6 ' 2 > 

with mean redshift (z) = zqF{A/ (3) / T{2>/ (3). We use the redshift distribution (6-2) 
with [3 = 1 and (z) = 3zq = 0.7 in our analysis. In Fig. 4 we show this redshift 
distribution p z (z) for f3 = 1 and (z) = 3zq = 0.7 and the relative lensing strength 
w(z) defined by Eq. (4-37) for the cluster redshift z d = 0.375. 

6.5. Results 

Figure 5 displays the two-dimensional /too-field for Abell 370 obtained using our 
galaxy sample and the redshift distribution (6-2) with (3 = 1 and (z) = 0.7. In Fig. 
6 we show the corresponding contour map of the /^-reconstruction. 

In the reconstruction, we have determined the two constants k<x>,a and Koo,b in 
the following way: (i) We employ the well-constrained mass model constructed by 
Kneib et al. 10 ) using the GLA and multiple images. The mean convergence within 
field A is calculated to be 0.832 using this model, (ii) Adopting this value for k^^, 

we determine k<x>,b such that the averages of Koo,a(%) arid Koo,b(%) within the 
overlapping region of the fields A and B are the same. The mean convergence in 
this overlapping region is 0.786, and the mean convergence in field B is 0.764. (iii) 
Joining together the two independent Koo-reconstructions along the diagonal of the 
overlapping quadrate, we obtain the resulting convergence map Koo(%), as was done 
in Ref 25). It can be seen in Figs. 5 and 6 that the discontinuity across this diagonal 
in our reconstruction is remarkably small. From the resulting K^-map, the mean 
convergence Roo inside the observed field of 0.65/i^~ 2 Mpc 2 is calculated to be 0.811, 
yielding a total mass in this field of M ~ 8.0 x 1O 14 /i^ o 1 M . 

The main features in the resulting mass map are summarized as follows: 

(1) Our mass reconstruction exhibits a clear bimodal feature in the central region 
where the two cD galaxies are located. The locations of the two mass maxima 
(^oo = 1-14) in the mass map coincide well with those of the two cD galaxies: The 
location of the northern mass clump is in good agreement with that of the northern 
cD galaxy, while the southern mass clump is located slightly offset from the southern 
cD galaxy, but is consistent with the location of the southern cD galaxy within a 
smoothing scale and a grid separation of ~ O.'l (37h^Q kpc). 

(2) The reconstructed mass map reveals a likely substructure (marked with * in Fig. 
6) located about 0.'3 (110 kpc) to the left relative to the center of the two cD galaxies. 

(3) An increase of Koo toward the lower-left corner in the lower-left quadrate is visi- 
ble. 

(4) The cluster-mass distribution is super-critical especially in the innermost region 
where the two cD galaxies and hence the two mass clumps lie. To be more precise, 
this distribution is super-critical for sources at an infinite or a sufficiently high red- 
shift. However, the reconstructed mass distribution is sub-critical for the GLA at a 
redshift of 0.724; w(z) = 0.45 for z = 0.724 (see Eq. (4-37)). 

Finally, it should be mentioned that these main features are insensitive to the as- 
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sumed redshift distribution p z (z) of field galaxies and the smoothing procedure which 
is needed in calculating the smoothed image ellipticity e. 




Fig. 5. The Koo-distribution for the cluster 
Abell 370 (z d = 0.375) reconstructed from 
the shear field shown in Fig. 3. The side 
length is about 2.'5 corresponding to 0.93 
Mpc at the cluster redshift z d = 0.375 for 
Ho = 50 km s _1 Mpc -1 and an Einstein- 
de Sitter universe. In the reconstruction we 
used the redshift distribution p z {z) defined 
by Eq. (6-2) with (3 = 1 and (z) = 0.7. Our 
reconstruction exhibits a clear bimodal fea- 
ture in the central region of the cluster. A 
likely substructure is visible near to the two 
cD galaxies. 
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Fig. 6. Contour plot of the resulting 

distribution shown in Fig. 5. The positions 
of the two cD galaxies are marked with +. 
The position of the northern mass clump 
agrees well with that of the northern cD 
galaxy. The southern mass clump is lo- 
cated slightly offset from the southern cD 
galaxy. The position of the extra substruc- 
ture is marked with *. The contours are 
stepped in units of 0.05. 



6.6. Discussion 

Focusing on the morphology of the cluster, we have investigated the mass dis- 
tribution of Abell 370 on the basis of weak lensing. We reconstructed the two- 
dimensional mass distribution of Abell 370 from the ellipticities of the faint galaxy 
images obtained using the HST/WFPC2 data. In the reconstruction, we used a 
total of 180 galaxy images consisting of 176 arclets, a radial arc and a GLA as a 
three-image system. Taking account of the small field-of-view of the data field, the 
strong-lensing features in the cluster and the fairly high redshift (z^ = 0.375) of the 
cluster, we applied the generalized mass-reconstruction scheme described in Section 
4.4. From the shear field of galaxy images, the mass distribution of the cluster was 
determined up to a one-parameter family of the global transformations which leaves 
the observed image distortions invariant. We employed the well-constrained mass 
model by Kneib et al. 10 ) based on strong lensing in order to infer the mass in the 
central region, from which we broke the mass degeneracy. Assuming a mean redshift 
of (z) = 0.7 for the field galaxies with i?675W £ (23,25), we have estimated the 
cluster mass within our data field of 0.65/^ 2 Mpc 2 to be ~ 8 x W u h^M Q . 

The resulting mass map also provides valuable information regarding the dis- 
tribution of dark matter. Our direct mass reconstruction exhibits a clear bimodal 



30 



K. Umetsu, M. Tada and T. Futamase 



feature associated with the two cD galaxies. Furthermore, our mass reconstruction 
reveals some other features: An extra substructure in the vicinity of the two cD 
galaxies and a mass condensation toward the lower-left boundary are visible. How- 
ever, since the accuracy of mass reconstructions tends to be worse near the boundary 
of the data field, the mass increase toward the boundary corner might be due to the 
systematic boundary effect and/or the fact that fewer galaxy images are used at 
the boundary corner in averaging image ellipticities. In addition, since the galaxy- 
number density in the upper-right WFC is quite small (see Fig. 2), the reconstructed 
mass distribution may be different from the actual one in that region. The result- 
ing mass distribution of the cluster is found to be super-critical for sources at a 
sufficiently-high redshift but sub-critical for the GLA at a redshift of 0.724. This 
result can be attributed to the smoothing procedure which is needed to reduce the 
noise due to the intrinsic ellipticities of background galaxies: Although the global 
mass distribution of a cluster can be obtained by locally averaging the galaxy-image 
ellipticities, such averaging will smooth out the galaxy-scale structure in the cluster. 
As a result, the reconstructed mass distribution will tend to be flatter than the orig- 
inal one. Hence, strong-lensing effects such as GLAs and multiple images should be 
combined with weak-lensing analysis in order to reproduce the mass distribution in 
the innermost region of the cluster. Alternatively, the galaxy-galaxy lensing analysis 
can probe the mass distributions in clusters on galaxy scales. 61 ) " 64 ) 

Finally, we should comment on our treatment of arcs: a GLA and a radial arc. 
In general, strong-lensing features such as GLAs and radial arcs strongly constrain 
the mass distribution of the deflector, especially in its central region where they are 
located. In fact, the reconstruction without taking into account the arcs exhibits a 
bimodal feature, but the density of the southern peak turns out to be lower than that 
of the northern peak. Thus the comparison between the resulting mass maps with 
and without the arcs reveals a remarkable influence of the arcs on the reconstructed 
mass distribution. 

§7. Summary 

Weakly- and coherently-distorted images of background galaxies induced by in- 
tervening galaxy clusters can be used to determine the mass distributions in the 
clusters on scales ~ 1 Mpc. Kaiser and Squires have derived an exact expression for 
the projected- mass distribution of the lensing cluster in terms of the gravitational 
shear which can be measured from the observed galaxy images. The original Kaiser 
and Squires technique was then generalized to include the strong-lensing regime and 
to optimize real observations. 

To this time, these direct mass-reconstruction methods based on shear analysis 
have been applied to a number of clusters. Now, the weak-distortion fields of galaxy 
images have become one of the most promising tools to probe the mass distribution of 
clusters. Moreover, the weak-shear field also provides a clue to the distant population 
of faint background galaxies: The weak-lensing analysis for a high-redshift cluster 
combined with a strongly-constrained mass model can place strong constraints on 
the redshift distribution of a distant population since the lensing strength depends 
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quite strongly on the source redshift through the distance ratio D^ s /D s (or w(z)) 
for a high-redshift cluster. The detection of significant lensing signals for some high- 
redshift clusters (z^ ~ 0.8) indicates that a substantial part of the faint galaxy- 
population must lie at sufficiently high redshifts (z > 1). 

In this paper we have applied a direct mass-reconstruction method to the cluster 
Abell 370 at a redshift of z d = 0.375 observed with the HST/WFPC2. Despite the 
small field-of-view and the irregular shape of the data field, our direct mass recon- 
struction of Abell 370 demonstrates the feasibility of weak-lensing analysis based on 
the shear field. However, there remains an uncertainty concerning the mass normal- 
ization of the cluster. Such an uncertainty occurs if the mass distribution of the 
cluster extends beyond the observed data field, in which case we cannot determine 
the cluster-mass distribution directly from the observed image distortions. We ex- 
pect that such a difficulty will be overcome by 8-10 m class telescopes with a wide 
field-of-view, such as the 8.3 m Subaru telescope with Suprime-Cam. 
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